Separability criterion for bipartite states and its generalization to multipartite systems 
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Quantum entanglement, introduced in the early days of quantum theory not only plays a vital role in difFerenti- 

■ ating quantum mechanics from classical mechanics, but also acts as key resource in quantum information processing. 

■ Thus the verification of quantum entanglement is fundamentally important in quantum information science. Since the 
mathematical description on quantum entanglement was first introduced by Werner in 1989 0, many entanglement 
measures have been proposed, such as distillable entanglement, entanglement cost Q, relative entropy of entanglement 

Schmidt number concurrence 0,01, negativity and so on @. Recently, the verification of multipartite 

■ entanglement and relevant issues are attracting more and more attention [lol fllj owing to the booming interest on 
quantum many-body systems [l2| . 

(Jj ■ The convexity requirement for quantum entanglement, which is in accord with the fact local operations with 
§ : classical communication (LOCC) c^n not create or increase quantum entanglement i , makes it very hard to verify 
the entanglement of mixed quantum states. Historically, positive partial transpose (PPT) criterion, proposed by Peres 
in 1996 Il3l . provides a necessary and sufficient condition for verifying the separability of 2 ig) 2 and 2 (g) 3 bipartite 
systems |14| . But for higher dimensional systems, how to distinguish entangled states from separable ones is still a 
challenging question. In this Letter, we introduce a group of symmetric operators to answer this question. It is shown 
~* \ that all these symmetric operators and their arbitrary linear combinations are entanglement witnesses. We use these 

■ symmetric operators to generalize Wootters' concurrence 0, 0] to high dimensional bipartite [isl Il6j and multipartite 
systems [l7| . More importantly, based on these symmetric operators the separability criterion for bipartite and 

. multipartite states is presented in arbitrary dimensions, which is a generalization of Wootters' separability criterion 

■ from the two-qubit case to a general case. 

, In the simplest two-qubit system, a pure state jV' >= '^ijNij2 >, where the subscripts "1" and "2" in the 

I ■ state vectors denote the two subsystems, is separable when it is a tensor product of a pure state in subsystem "1" and 
psj '. a pure state in subsystem "2" , thus the four coefficients £ {1, 2}) have to satisfy the condition 011022 = 112021. 

^-H ■ In Ref. , Wootters introduced a "spin-flip" operator dy (Pauli matrix) in the standard basis to verify the separability 
of a two-qubit pure state. That is, a two-qubit pure state \ip > is separable iff < ^*\ay® ay\ilj >= 0. In the following, 
we call the state vector \ip > and its complex conjugate > orthogonal with respect to the operator O, or simply 
O-orthogonal, if they satisfy the relation < 0. Wootters' separability criterion for a two-qubit pure state 

can then be stated as, it should be O-orthogonal to its complex conjugate with the operator O = ay ay. Here 
the antisymmetry of the "spin- flip" operator ay, i.e. crj — —ay (the superscript "T" stands for matrix transpose 
hereafter), plays the essential role in verifying the separability of a two-qubit state. This O-orthogonal relation 
between a separable two-qubit pure state and its complex conjugate is violated by all entangled two-qubit states. 

Now we consider a general bipartite case, with the two subsystems having dimensions Di and D2, respectively. 
If a bipartite pure state \^ >= X^iii "ij > is separable, the four coefficients a^j, a^j/, at'j, and Qi'j' with 

i, i' € {1, 2, • • • , Di} and j, j' £ {1, 2, . . . , D2}, have to satisfy the condition aijai'j' = aij'ai'j, so that this bipartite 
pure state can be formulated as a tensor product. This condition, similar to the above two-qubit case, is equivalent 
to the following O^*'* ^-orthogonal relation, 

< >^0, (la) 

with 

0(M':j-./) ^ [^("') _ {a[''''Y] ® [ai'''"> - {ai"">f]. (lb) 

Here the operator a^'^ ■* (ctj"''"' ^) in the subsystem "1" ("2") has only one nonzero element "1" located at row i (j) and 
column i' {j') in the standard basis, no matter what dimensions it has. The two antisymmetric (or skew-symmetric) 
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operators, [trj*'' ■* — (ci''* ■')"^] and [(T2"''"' ■* ~ ('^2"''"' which play the same role as the antisymmetric operator ay in 
the two-qubit case, establish a relation among the four coefficients aij, aij' , aiij, and atiji. 

Since the superscript numbers € {1,2,- •• ,-Di} and j,j' G {1, 2, . . . , D2} in the above Eq.([T]) can be chosen 
arbitrarily, the O^*'* ■■''^ ■'-orthogonal relation ([Ta]) for a separable bipartite pure state, holds valid for all possible 
operators O^*'* '^'^ ■* in the form (|lbp . Owing to the linear dependence of the relation ([Ta]) on the operators O'*'' ■■''^ \ 
we even can replace the operators cr^*'* ^ and (T2'^'"' by a random matrix u^"^™'' in subsystem "1" and a random matrix 
^(rdm) subsystem "2" to construct, 

Q(Semi-rdm) ^ [^(rdm) _ (^("-drnj^Tj ^ [^(rdm) _ (^('■dm))Tj^ (2) 

and the 0(Semi-rdm)_Qj.^j^QgQj^g^l relation, < V* 10^^''™"'''^™^ I?/- >= 0, stiU holds true if only the pure state \tp > is 
separable. Thus we can conclude that a bipartite pure state in arbitrary dimensions is separable iff it and its complex 
conjugate are orthogonal with respect to all symmetric operators formulated as a tensor product of two antisymmetric 
operators in the two subsystems. To verify the separability of a given bipartite pure state, we only need to pick out 
a finite number of linearly independent symmetric operators and check the above 0(^°™'"'''^™^-orthogonal relation. 
A violation of the 0*^^''™'"''^™-'-orthogonal relation by any operator (^(Scmi-rdm) ^j^Q^g]^ declare the existence of 
quantum entanglement. 

Now we use a complete set of linearly independent symmetric operators O*^*'' '^'^ \ which are directly constructed 
by ^Di{Di — 1) linearly independent antisymmetric operators [a[^'^ •* — (crj*'' •')^] in subsystem "1" and \D2{D2 — 1) 

linearly independent antisymmetric operators ^ ^{(^2 '^ ^Y"] subsystem "2" , to define the entanglement measure 
of a bipartite pure state IV' >, 

E e' E i<^i[ar'-(ar^n®[4"''-(4""''n^ o) 

2=1 i'=i+l j = l j'=i+l 

where the superscript "(2)" denotes the bipartite system. This operator-based entanglement measure for bipartite 
pure states, which we call concurrence in the following, turns back to Wootters' concurrence for two-qubit systems, 
and satisfies the basic requirements for a good entanglement measure. For example, (i) it presents zero result for 
all tensor product states, and positive results for entangled states; (ii) it remains invariant under all local unitary 
transformations, and so on. In fact, this O concurrence is the same as measured in the following way, 

C(2)(|^>) = y'2-tr(p?)-tr(pi), (4) 

with pi {P2) being the reduced density matrix of subsystem "1" ("2"), and "tr()" denoting matrix trace. Since the 
equality tr(pf ) = tr(p2) holds for bipartite pure states, our entanglement measure is equivalent to the entanglement 
measure proposed in Ref . [l6| . which is called / concurrence owing to its connection to the universal inverter. It is 
very interesting that two measures present the same result for the bipartite entanglement, though they are defined in 
two different ways and based on different ideas. 

This O concurrence can also be extended by convex roof to measure the entanglement of mixed bipartite states. 
At present, we only concern the separability condition of mixed bipartite states, rather than the evaluation of their 
entanglement. According to Werner's criterion a mixed bipartite state p is separable only when it can be written 
as a mixture of separable pure states, 

p = Y, p^m X V'i I ® W2 xi^ii (5) 

i-J 

with the probabilities Pij > and Pij ^ 1- However, every mixed quantum state has countless types of pure- 
state decomposition. Given a density matrix, it is usually very difficult to prove the existence or nonexistence of a 
decomposition composed of product states. In other words, we can easily construct a separable mixed quantum state 
by mixing a group of separable states in the above way ([S]). However, it is very difficult to recover this group of 
separable states, provided only the density matrix. 

We begin our investigation on the separability of a mixed Di <Si D2 state with a Hermitian matrix y^, which is 
the square root of the density matrix p. Every pure-state decomposition of the density matrix p is connected to a 
D12 X M {D12 = Di X D2 < M) right-unitary transformation U through the relation, W = . Each column 
vector of the matrix W , the ith column vector denoted as \Wi > hereafter, is a "subnormalized" pure state in the 
[/-decomposition, and the original density matrix can then be written as p = X)i=i Now we suppose 



c<"'(W >) 
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that the density matrix p represents a separable mixed state, which means there exists at least one U -decomposition 
so that all pure states in this decomposition are tensor product states. That is to say, every column vector \Wi > is 
O^^'^"""''^'"-' -orthogonal to its complex conjugate, with the operator (3(SGmi-rdm) (^ggng^j Eq.(l2|). 

Since all column vectors \W, > satisfy the 0'-^°"''-"^"'^ -orthogonal relation, < W* \0''^'"^"'^"'^W, >= 0, for a 
particular [/-decomposition of the separable mixed state p, the symmetric matrix M^^0(^°™'"'''^™)M^ is a hollow matrix 
whose diagonal elements are all equal to zero. Equivalently speaking, a mixed bipartite state p is separable, iff there 
exists a right-unitary transformation U, so that the symmetric matrix SU , with 

is a hollow matrix for all operators (^(Somi-idm) ^^^^ form ©. Here we emphasize the above symmetric matrix S 
depends only on the given density matrix and an operator constructed by two random matrices. Now a new question 
arises, given a symmetric matrix 5, whether there exists such a right-unitary transformation U so that U'^ SU is a 
hollow matrix? 

We can find out the answer by using Woottcrs' method in Ref.Q. A symmetric matrix S can be transformed to a 
hollow matrix in the way of U'^ SU, only when its maximal singular value is no larger than the sum of the rest singular 
values. Supposing {A^} are singular values of the symmetric matrix S in decreasing order, the above statement is 
equivalent to the following condition, 

Ai<^A,. (7) 

1=2 

On the contrary, if the above condition ([7]) is not satisfied for even one particular operator in the form ([2]), e.g. 
O(particuiar) ^ J^f^i^ T,?L^+lIlf=l^ T,f'lj+i c^.^' -.j .j' O^''' the mixcd statc p must be entangled. Please see 

Ea. (jlbp for the definition of the operator O^*'* '^'^ \ Without loss of generality, we assume the first coefficients Ci2:i2 
has the maximal modulus among all the complex coefficients {ci^i'-jj/}. In this case, the average O concurrence of 

the pure states, i.e. column vectors {\Wk >}, in an arbitrary [/-decomposition, C^\p) = X^fcli C'^^Hl^fe is 

larger than zero, because C'lf\p) > . ^ (Ai — J^fll -^j)- given a mixed bipartite state p, if the condition ([7]) 

Y|ci2:12| -' 

is violated by even one operator in the form ([2]), this mixed state is an entangled state. One of main conclusions 
in this Letters can now be concluded, a mixed bipartite state in arbitrary dimensions is separable only when the 
condition ([7]) for the singular values of the symmetric matrix ^ is generally true for all operators (^(Scmi-rdm) ^-^^ 
form (12) . This is a necessary and sufficient condition. On the one hand, if the condition ([7]) is valid for all operators 
Q(Semi-rdm) ^-j^g form there must exist a unitary transformation U independent on the operator (^(Somi-rdm) 

(because it is constructed by a random matrix a^'^"^^ in subsystem "1" , and a random matrix CTj'^'^™^ in subsystem 
"2" ), to satisfy all 0^^°™'"'''^™^-orthogonal relations in the corresponding decomposition, and such a bipartite state 
is separable; On the other hand, any violation of this condition, just as we already shown, is enough to declare the 
existence of quantum entanglement. In other words, every operator (^(Somi-idm) ^^i^ form ([2]) is a witness [l^ of 
bipartite entanglement under the violation of the condition Q. 

The singular values {A;} of the symmetric matrix S ([6|) are equal to the square root of the eigenvalues of the matrix 
S^S [ll], and as well the square root of the eigenvalues of the following matrix, 

p ^(Scmi-rdm) 

which is very similar to the form used in the concurrence paper , and can be considered as its generalization in high 
dimensional systems. 

At present, we can not verify the separability of a bipartite state in high dimensions through finite tests on the 
condition ([7]), unless it is a pure state. But we here provided an efficient way for verifying the entanglement of a 
bipartite statc, because any operator (^(Scmi-rdm) ^j^g form ^ is an entanglement witness under the violation of 
the condition ([T]). Our numerical results show that the larger entanglement a bipartite state contains, the fewer tests 
it usually requires to find a violation by randomly generating the matrices (tJ"^™'' and ctj"^™'' in the form ([2]). 

Now we come to the multipartite case. A multipartite pure state is said to be entangled only when it cannot 
be written as a tensor product of the states in the subsystems. This definition does not distinguish between "truly 
multipartite" entanglement and low-partite entanglement (20| . For example, both the Greenberger-Horne-Zeilinger 
state and the W state in the three-qubit system are tripartite entangled states according to the above definition, 
but only the former one has nonzero 3-tangle for "truly tripartite" entanglement jTo| . Our following discussion on 
the separability criterion of multipartite states is based on this definition of entanglement, no matter what kind of 
entanglement, "truly multipartite" entanglement or low-partite entanglement, is contained. 



^*Q(Scmi-rdm)^ (8) 
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For a separable multipartite pure state, 

\ip >^ It/Ji >(E)\tp2 ><E)- ■ ■ (E}\ipN >, (9) 

each subsystem is separable with the rest subsystems as a whole. The inverse situation is also true. Among a 
multipartite system, if every subsystem is separable with the rest subsystems as a whole, this multipartite pure 
state is separable. That is to say, multipartite entanglement can be featured by bipartite entanglement. Thus the 



entanglement degree of a multipartite pure state can be quantified by C^^'^{\Tp >) ~ y 5 ^k=i I^^^HIV'^fe >)Pj where 
the superscript "(A^)" denotes A^-partite system, C*^^^ is the bipartite entanglement defined in Eg.®, and > is 
the bipartite version of the quantum state {ip > in Eq.® by considering (N — 1) subsystems, excluding the fcth one, 
as a whole. Although some other bipartite entanglement, e.g. C**-^-* (lV'(i2)(T2) >)' where the first part (12) is composed 
of the two subsystems "1" and "2", and the other part (12) is composed of the rest (TV — 2) subsystems, can also 
be included to define other types of multipartite entanglement, for example in Ref. |2l| . the participation of these 
terms only affects the evaluation of entanglement, but does not change the separability criterion. It means we can 
measure the multipartite entanglement in a simpler way, where every involved bipartite entanglement is associated 
with one single subsystem and the rest as a whole. The above entanglement measure for multipartite states can also 
be described in terms of operators. 



c<">(|* » = 



\ 



^EE---E i<^ior-"^'i'^>^ (10) 



where the symmetric operator, 



is associated the bipartite system composed of the fcth subsystem and all the other subsystems as a whole, Dk is 
the dimension degree of the fcth subsystem and T.c. means matrix transpose. This is generalized O concurrence 
for multipartite pure states, which, inherited from bipartite entanglement, remains invariant under local unitary 
transformations. Under this measure, the maximally-entangled D-dimensional A^-partitc Greenbergcr-Horne-Zeilinger 

state |Gi/Z(^) >= X^i^i Ni*2 • • • > has entanglement ^ N{1 — -p), and the generalized A^-qubit W state 

>= -^(|100 • • • > +|010 • • • > H h |00 • • • 1 >) has entanglement ^2{1 - jj). Similar to the bipartite case, 

the entanglement measure (|10|) has an equivalent but much simpler version. 



C('"(|*>) 



N 



N-J2tr{pl), (12) 



fe=i 



where Pk stands for the reduced density matrix of the fcth subsystem. This measure is a generalization of the / 
concurrence [l6j for multipartite systems. 

Since the present entanglement measure for multipartite systems is based on the bipartite entanglement between 
single subsystems and the rest as a whole, we can directly generalize some conclusions from bipartite case to the 
present multipartite case. For example, a multipartite pure state is separable only when it is orthogonal to its 
complex conjugate with respect to all operators in the form (jll[) . and also their extended version, 

^(Semi-rdm) ^ ^(rdm) ^ , , , ^ j^(rdm) _ (^("'drnj^T] ^ . . . ^ ^(rdm) ^ r^^^_^ ^^^^ 

where cr^"^™-* is a random Dk x Dk matrix in the fcth subsystem. 

For a separable multipartite state p, there exists a decomposition, W = ^/pU, where each "subnormalized" pure 

state \Wi > (ith column vector of W) is a tensor product state in the form (|9]), thus is 0^^''™'''^^™''-orthogonal to its 

complex conjugate. Similar to the above case, the matrix W'^ O^^'^^'^'^^'W has to be a hollow matrix, which means 
the singular values {A^} in decreasing order of the symmetric matrix, 

(14) 

have to satisfy the condition ([7|). Now we conclude the separability criterion for a mixed multipartite state, that is, 
the condition ([7]) for the singular values of the above matrix p4l) holds valid for all symmetric operators in the form 
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(fT3|) . A violation by any operator in the form (|13p is enough to declare the existence of multipartite entanglement. 
So every operator jj^ ^j^^ form (|13p is a witness of multipartite entanglement under the violation of the 

condition ([7]). This huge class of entanglement witnesses can help us efficiently verify multipartite entanglement, even 
for mixed multipartite states. Here we also note that the singular values {A^} of the symmetric matrix (|14p are the 
square root of the eigenvalues of 5^ 5, and as well the square root of the eigenvalues of the following matrix, 

O' (Scmi-rdm) 
■ I kk 

which is directly determined by the density matrix p and an operator o'"^"^^ '"^"^^ constructed by random matrices in 
the way p3|) . 

To summarize, a group of symmetric operators with only two nonzero matrix elements are introduced to carry out 
the separability criterion of a general quantum state. Every symmetric operator, formulated as a tensor product of an 
antisymmetric operator in one subsystem and an antisymmetric operator in other subsystems as a whole, and their 
arbitrary linear combinations are entanglement witnesses for multipartite quantum states. Based on these symmetric 
operators, Wootters' separability criterion for two-qubit states is generalized to multipartite systems in arbitrary 
dimensions. How to distinguish all mixed entangled states from separable ones within a finite number of tests is the 
next challenging question worthy of investigation. Its answer might directly bring us a computable entanglement 
measure for mixed multipartite states. 

This work was supported by the national Natural Science Foundation of China under Grant Nos. 11174118, 11174026 
and 11264018, and the Natural Science Foundation of Jiangxi Province, China under Grant No. 20114BAB212003. 
We also thank Prof. J. P. Dowling for helpful discussions. 

Note added.- After completion of this work, we became aware of a related work by Chen, Ma, Giihne, and Sevcrini[22l|. 
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